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Abstract — In  this  study,  the  Kantorovich-Vlasov  method 
has  been  applied  to  the  flexural  analysis  of  simply 
supported  Kirchhoff  plates  under  transverse  uniformly 
distributed  load  on  the  entire  plate  domain.  Vlasov  method 
was  used  to  construct  the  coordinate  functions  in  the  x 
direction  and  the  Kantorovich  method  was  used  to  consider 
the  assumed  displacement  field  over  the  plate.  The  total 
potential  energy  functional  and  the  corresponding  Euler- 
Lagrange  equations  were  obtained.  This  was  solved  subject 
to  the  boundary  conditions  to  obtain  the  displacement  field 
over  the  plate.  Bending  moments  were  then  obtained  using 
the  moment  curvature  equations.  The  solutions  obtained 
were  rapidly  convergent  series  for  deflection,  and  bending 
moments.  Maximum  deflection  and  maximum  bending 
moments  occurred  at  the  center  and  were  also  obtained  as 
rapidly  convergent  series.  The  series  were  computed  for 
varying  plate  aspect  ratios.  The  results  were  identical  with 
Levy-Nadai  solutions  for  the  same  problem. 

Keywords — Kantorovich-Vlasov  method,  Levy-Nadai 
method,  convergent  series,  total  potential  energy 
functional,  Euler-Lagrange  differential  equation. 

I.  INTRODUCTION 

Plates  are  widely  used  in  many  engineering  structures  like 
roof  and  floor  slabs,  bridge  deck  slabs,  foundation  footings, 
bulkheads,  watertanks,  turbine  disks,  spacecraft  panels  and 
ship  hulls.  Plates  may  be  subjected  to  static  flexural  loads, 
dynamic  loads  or  inplane  loads,  resulting  in  three  types  of 
analysis,  namely  - static  flexural  analysis,  dynamic  flexural 
analysis  and  stability  analysis  [1],  The  deformation  of  plates 
is  usually  defined  as  the  deformation  of  the  middle  surface, 
which  is  a plane  surface  equidistant  from  the  top  and 
bottom  faces  of  the  plate  [2],  The  thickness  of  the  plate  h,  is 
the  distance  between  the  top  and  bottom  faces,  measured  in 
a direction  normal  to  the  middle  plane. 

Generally,  plates  are  classified  into  three  groups:  thin 
plates,  moderately  thick  plates  and  thick  plates,  depending 
upon  the  thickness  of  the  plate  [2],  Plates  are  also  classified 


depending  on  the  nature  of  deformation  and  material 
properties  as  plates  with  small  deformations,  plates  with 
large  deformations,  orthotropic  plates,  isotropic  plates  and 
anisotropic  plates.  They  can  also  be  classified  according  to 
their  shapes  as  rectangular,  circular,  triangular,  skew  and 
elliptical  plates. 

Thin  plates  are  plates  for  which  the  thickness  is  very  small 
when  compared  with  the  smaller  of  the  in  plane  dimensions. 
For  most  practical  applications,  thin  plates  have  a ratio  of 
thickness  to  the  smaller  in  plane  dimension  less  than  0.05 
[2],  The  plate  deformation  is  considered  small  if  the 
maximum  deformation  is  less  than  or  equal  to  its  thickness. 
The  most  widely  used  plate  theory  is  the  Kirchhoff-Love 
plate  theory,  also  called  the  classical  thin  plate  theory  [3], 
The  Kirchhoff-Love  plate  theory  is  based  on  the  Kirchhoff- 
Love  hypothesis  which  makes  assumptions  similar  to  those 
in  the  Bernoulli-Navier  hypothesis  used  in  the  theory  of  thin 
(shallow)  beams. 

For  large  deflections,  the  flexure  of  the  plate  is 
accompanied  with  inplane  stresses  which  depend  not  solely 
on  the  deflection  but  also  on  the  type  of  edge  supports  [2], 
The  incorporation  of  inplane  stresses  in  the  behaviour  of 
plates  results  in  non-linear  differential  equations  which 
present  theoretically  rigorous  demands  for  solutions. 

For  a good  number  of  engineering  applications,  the 
Kirchoff-Love  theory  gives  sufficiently  accurate  results,  the 
accuracy  however  reduces  with  increased  plate  thickness, 
with  localized  loads,  and  in  problems  of  stress 
concentrations  around  openings  [2],  [4],  Again,  the 

Kirchhoff-Love  plate  theory  sometimes  uses  an 
approximate  incomplete  set  of  boundary  conditions  which 
results  in  errors,  especially  near  the  boundary  regions  as 
well  as  in  reactions  [2],  [4], 

Thick  plate  problems  can  only  be  formulated  and  solved 
using  the  three  dimensional  theory  of  elasticity.  For  many 
engineering  problems,  solutions  based  on  the  three 
dimensional  theory  of  elasticity  present  mathematically 
rigorous  demands  for  solution,  and  thus  far,  only  simple 
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plate  problems  have  been  solved  using  this  theory  [2],  [4], 
[5],  [6],  [7], 

Several  other  theories  have  been  proposed  to  describe  the 
behaviour  of  plates  under  loads.  They  include  - Reissner’s 
sixth  order  plate  theory  [8],  [9]  for  moderately  thin  plates; 
Von  Karman  plate  theory  for  plates  under  large 
deformation;  Mindlin  plate  theory  for  moderately  thick 
plates,  shear  deformation  plate  theories  [10,  11,  12], 

The  focus  of  this  paper  is  the  Kirchhoff  plate  theory. 
Several  methods  have  been  used  to  analyse  plates.  They  are 
broadly  classified  as  analytical  methods  and  numerical 
methods.  Analytical  methods  include  Navier  and  Levy 
methods  [13].  Numerical  methods  which  seek  approximate 
solutions  to  the  plate  problem  include:  variational  methods, 
finite  element  and  finite  difference  methods  and  finite  strip 
methods  [14,  15,  16]. 

Kirchhoff-Love  hypothesis 

The  Kirchhoff-Love  hypothesis  for  the  small  deflection 
theory  of  thin  homogeneous  elastic  plates  are  [2]: 

i.  the  deflection  of  the  middle  surface  is  small 
compared  to  the  thickness 

ii.  the  middle  surface  is  unstrained  subsequent  to 
bending,  and  hence  is  a neutral  surface 

iii.  plane  sections  initially  normal  to  the  middle  surface 
remain  plane  and  normal  to  the  middle  surface  even 
after  deformation 

iv.  the  stress  (cl)  normal  to  the  middle  surface  is  very 
small  compared  to  the  other  two  normal  stress 
components  ( Gxx  and  G ) and  hence  G_  can  be 
neglected. 

These  assumptions  reduce  the  three  dimensional  elasticity 
problems  to  two  dimensions. 

II.  RESEARCH  AIM  AND  OBJECTIVES 

The  aim  of  this  study  is  to  apply  the  Kantorovich-Vlasov 
method  to  the  analysis  of  simply  supported  thin  rectangular 
plates  under  uniformly  distributed  transverse  loads.  The 
specific  objectives  are: 

i.  to  use  the  Vlasov  procedure  to  derive  suitable 
coordinate  shape  functions  for  the  simply  supported 
plate  in  the  x coordinate  direction. 

ii.  to  obtain  the  total  potential  energy  functional  Id  for 
the  plate  bending  problem  based  on  the  Vlasov  shape 
functions  derived. 

iii.  to  obtain  the  Euler  Lagrange  differential  equations 
for  the  total  potential  energy  functional  derived. 


iv. 


v. 


VI. 


to  solve  the  Euler-Lagrange  differential  equations 
subject  to  the  boundary  conditions  of  the  plate  in  the 
y coordinate  direction. 

to  use  the  bending  moment  displacement  relations 
and  obtain  the  bending  moment  expression, 
to  obtain  the  maximum  values  of  deflection  and 
bending  moments  at  the  centre  of  the  plate  for 
various  aspect  ratios. 


III.  THEORETICAL  FRAMEWORK 

The  strain  energy  of  a plate  for  linear  elastic  behaviour  is 
given  by  the  volume  integral: 

hi  2 


U = \\\  1 atiev  dzdxdy 


R-h/2 

0) 

where  G(/  are  stresses,  S;y  are  strains,  h is  the  plate 

thickness,  R is  the  two  dimensional  domain  of  the  plate  on 
the  xy  plane. 

From  the  stress-strain  relations,  the  strain  energy  expression 
can  be  expressed  in  terms  of  strains  and  stresses  as  follows: 


/i/2 


u = 


2(1  E ) R-h/2 

(2) 


J J J [eL  + 2fl£xc£xv  + syy  + 2(!  _ 


u = — 

2 E 


l-  JJJ  {(°x  + + a2z  - 2\i(axay  + axaz  + ayo 


+ 2(1  - |4)(T^  + Txz  + t y jjdxdydz 
(3) 

where  Sxr,  Syy  are  normal  strains,  (Jxx,  Cfvy  are  normal 
stresses,  £ fv  is  the  shear  strain,  O vv  is  the  shear  stress,  ft 

is  the  Poisson’s  ratio,  E is  the  Young’s  modulus  of 
elasticity,  and  V is  the  three  dimensional  region  of  the  plate. 
If  the  distributed  transverse  loads  p(x , y)  are  assumed  to  act 
on  the  middle  surface  of  the  plate  and  w(x,  y)  is  the  vertical 
displacement  function  of  the  plate  middle  surface,  the 
potential  energy  of  the  external  distributed  loads  is  given 
by: 
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Ve  = -jj  p(x,  y)w(x,  y)dxdy 


(4) 

For  small  deformations,  the  strain-displacement  relations 
are: 


d2w 


^ XX  Z- 


dx2 


(5) 


£yy  = -Z 


d2w 


8y2 


(6) 


£xy  =~Z 


d2w 

dxdy 


(7) 


where  w(x,  y)  is  the  transverse  deflection. 
Using  the  stress-strain  laws, 

n = u+v. 


°xx  - - 


Ez 


(8) 


Ez 


(1-r) 


f X2  a2  'l 

o w o w 

9 f ft T 

v dx2  dy 2 j 


( 0,2  T 

o w o w 

— T + l11 — V 
8y2  8x2 


(9) 


x*y  . . 2 


Ez  d2w 


(1  — ju  ) dxdy 

(10) 

= ^zy  = Vzz  = 0 

(ID 

where  E is  the  Young’s  modulus  of  elasticity,  |i  is  the 
Poisson’s  ratio,  T„V,T_,  and  T„,  are  shear  stresses,  a„ 

Ay  ZZ 

is  normal  stress. 

The  total  potential  energy  functional  becomes: 
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(12) 


where  D = 


Elt 


12(1 -Hz) 

D is  the  bending  rigidity,  h is  the  plate’s  thickness. 


(13) 


IV.  APPLICATION  OF  THE  KANTOROVICH- 
VLASOV  METHOD  TO  THE  BENDING 
ANALYSIS  OF  SIMPLY  SUPPORTED 
KIRCHHOFF  PLATES 

Consider  the  rectangular  Kirchhoff-Love  plate  shown  in 
Figure  1. 


r i iiiiiii  . 

7W 

Fig.l:  Rectangular  Kirchhoff-Love  plate  under  uniformly 
distributed  load. 


The  coordinate  axis  is  chosen  as  shown  in  Figure  1 , in  order 
to  take  advantage  of  the  symmetry  of  the  problem  in  the  y 
coordinate  direction.  The  plate  is  subject  to  a transverse 
distributed  load  of  intensity  q over  the  whole  plate  domain 

which  is  given  by  0 < X < a,  — ^4  < y < “/ 


4.1  Vlasov  procedure  for  finding  the  displacement 
shape  function  in  the  x-direction 

Following  Vlasov  procedure,  the  displacement  shape 
function  in  the  x-direction  are  the  eigenfunctions  of  an 
Euler-Bernoulli  beam  on  the  x-axis.  The  differential 
equation  of  free  vibration  of  Euler-Bernoulli  beams  is 


E!  d4X(x,  t) 
dx 4 


+ in 


d2X(x,t) 
dt 2 


= 0 


(14) 
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where  El  = flexural  rigidity  of  beam.  111  = mass  per  unit  length  of  beam,  X(x,  t ) are  displacement  functions  in  the  x direction. 
For  harmonic  vibrations,  X will  vary  with  space,  x and  time,  t variables  as: 

X(x , t)  = f(x)  sin  to nt 

where  00;i  is  the  natural  frequency,  and  /(v)  is  the  displacement  shape  function  in  the  ^-direction. 

Then  the  differential  equation  of  free  vibration  becomes: 


/ivw-^/w 

El 


\ 


sin  cot  = 0 


y 


fiv(X)-t^f(x)  = 0 


where  f"  (x) 


II 

d4f(x) 

dxA 
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(15) 


(16) 


or 


r\x)-^f(x)= o 
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where 

El  a4 

The  general  solution  of  Equation  (17)  is: 


(3x 


. . px  Px  . . Px  . 

/ (x)  = q sin J h c2  cos  — + c3  sinh  — + c4  cosh 

a ~ a a a 


(18) 


where  ci,  ci , C3  and  C4  are  integration  constants,  which  can  be  found  from  the  boundary  conditions,  and  P is  the  root  of  the 
characteristic  equation. 

The  eigenfunctions  f for  the  mth  mode  of  vibration  are  obtained  as 


( 


/mU0  = 


clm  sin  -s—  + c2m  cos  + c3m  sinh  + c4m  cosh  -2L- 


v 


a j 


(19) 


a a a 

For  simple  supports  at  x = 0,  x = a the  boundary  conditions  are 

/m(0) = fm(°)  = 0 

/*(«)  = /»  = 0 

Using  the  boundary  conditions  in  Equation  (20),  the  following  algebraic  homogeneous  equation  is  obtained: 

0 
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clm 

-1 

0 
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(21) 


The  solution  yields 


fm(x)  = clm  sin 


mux 


a 


m = 1,  2,  3, 


Thus,  the  eigenfunction  is  obtained  as 

„ . . . mux 

fm(x)  = Sin 


(22) 


(23) 


a 


where  m = 1,  2,  3,  4,  ... 
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4.2  Kantorovich  procedure  to  obtain  the  total  potential  energy  functional 

Following  Kantorovich  procedure,  the  plate  deflection  is  considered  as 


w(x,y)=  2#OOsin 


m= 1 


mux 


a 


(24) 


where  g(y)  is  an  unknown  function  of  the  independent  coordinate  variable  v.  Substitution  of  Equation  (24)  into  the  total  potential 
energy  functional  yields  after  simplification: 
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The  integrand  in  I I is: 
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Or,  by  division  by  Da, 
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where  F,  and  Fj  are  integrands  of  II  which  depends  on  y,  g(y ) and  g"(y) 


(29) 


(30) 


4.3  Euler  Lagrange  differential  equation 

The  total  potential  energy  functional  is  minimized  when 

sn  = o 


This  corresponds  to  the  Euler  Lagrange  differential  equation: 


dF 
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From  Equation  (32),  the  Euler-Lagrange  differential  equation  is 


(31) 

(32) 
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g,v(y)-  2 
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(33) 


4.4  Solution  of  the  Euler-Lagrange  equation 

The  Euler-Lagrange  differential  equation  is  solved  using  the  method  of  trial  functions,  to  obtain  the  homogeneous  solution  as: 

. . , mny  „ mny  . , mny  _ mny  , mn y . , mny 

g(y ) = Am  cosh  — - + Bm  — - sinh + Cm cosh + Dm  sinh (34) 

a a a a a a 

where  Am,  Bm,  Cm  and  D„,  are  the  four  constants  of  integration. 

The  particular  solution  gp(y ) for  uniformly  distributed  load  q(y  ) = qa  is  given  as 

«PW  = 7^Trr)  <35> 

(mnY  D 

This  problem  is  symmetric,  hence  Cm  = Dm  = 0 since  the  terms  of  g(y)  associated  with  them  are  odd  functions  of  y.  Thus  the 
general  solution  is: 
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The  two  constants  of  integration  C2m  and  csm  are  found  using  the  geometric  boundary  conditions  W | X,  y = + j = 0 and  the 


force  boundary  conditions  y = + j = 0 


From  the  force  boundary  conditions 
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From  the  geometric  boundary  conditions. 
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Solving  simultaneously,  the  constants  are  obtained  as: 
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This  completely  determines  the  deflection. 


Center  deflection 

The  maximum  deflection  occurs  at  the  plate  center  (x  = all , y = 0),  and  is  given  by: 
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m = 1,  3,  5,  7,  ... 
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Bending  Moment  Expressions 

From  the  bending  moment  displacement  relations,  the  bending  moments  are  obtained  as: 
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Maximum  bending  moments  occur  at  the  plate  center  where  x = all,  y = 0 and  are  given  by: 
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c21  = -2. 742459 

Dn5 


C31  _ 


2 q0a4 


Dn5  cosh  — 


c31  = 0. 797074 

Dn5 

Mxx  = 0.051668<70a' 


(44) 


(45) 


(46) 

(47) 

(48) 

(49) 

(50) 

(51) 

(52) 
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For  m = 1,  3;  b/a  = 1 

Mxx  = 0.0472387g0a2  (53) 

For  m = 1,  3,  5;  b/a  = 1 

Mxx  = 0.0419q0a2  (54) 

Similarly  for  m =1, 

4 

w =0.00411^  (55) 

c D 

For  m = 1,3 

4 

w = 0.00406^-  (56) 

c D 

The  expressions  for  the  maximum  deflection  and  the  maximum  bending  moments  are  evaluated  at  the  plate  center  for  various 
values  of  the  aspect  ratio,  and  tabulated  in  Table  1. 


Table. 1:  Coefficients  for  maximum  deflection  and  maximum  bending  moments  in  simply  supported  plates 


b/ 

/ a 

4 

q0a 

Wmax  = a 

Mxx  = Px^0«2 

Myy  = Pvv^o"2 

1 

0.00406 

0.0477 

0.0479 

1.2 

0.00564 

0.0626 

0.0501 

1.4 

0.0753 

0.0506 

1.5 

0.00772 

0.0812 

0.0498 

1.6 

0.0862 

0.0493 

1.7 

0.00883 

0.0908 

0.0486 

1.8 

0.0948 

0.0479 

2 

0.01013 

0.1017 

0.0464 

3 

0.01223 

0.1189 

0.0404 

4 

0.01282 

0.1235 

0.0384 

V.  RESULTS  AND  DISCUSSIONS 

The  Kantorovich- Vlasov  method  has  been  successfully 
applied  to  the  bending  problems  of  simply  supported 
rectangular  Kirchhoff-Love  plates  subjected  to  uniformly 
distributed  transverse  loads  over  the  entire  plate  domain. 
The  problem  of  bending  of  Kirchhoff-Love  plates  under 
transverse  distributed  load  was  presented  in  variational 
form  and  the  total  potential  energy  functional  for  Kirchhoff- 
Love  plates  presented  as  Equation  (2).  Vlasov  procedure 
was  adopted  to  construct  the  displacement  coordinate 
functions  in  the  jt-di  recti  on  as  Equation  (23).  Kantorovich 
method  was  adopted  to  consider  the  displacement  function 
of  the  plate  as  the  series  given  by  Equation  (24).  This 
assumed  displacement  function.  Equation  (24)  was 
substituted  into  the  total  potential  energy  functional  to 
obtain  the  total  potential  energy  functional  expressed  by 
Equation  (28),  where  the  integrand  was  found  to  depend  on 
the  functions,  g(y).  and  g"(y).  The  Euler-Lagrange 


differential  equation  corresponding  to  the  extremization  of 
n was  obtained  as  Equation  (33),  a fourth  order  ordinary 
differential  equation  in  g(y ) and  its  derivatives.  The  Euler- 
Lagrange  differential  equation  was  then  solved  using  the 
trial  function  method  to  obtain  the  solution  as  Equation  (36) 
after  taking  consideration  of  symmetry  of  plate  and  loading. 
Geometric  (displacement)  and  force  boundary  conditions 
were  used  to  obtain  the  constants  of  integration  as 
Equations  (40)  and  (41).  The  maximum  deflection  was 
found  to  occur  at  the  plate  center,  in  line  with  symmetry  and 
the  maximum  deflection  was  obtained  as  Equation  (43). 
Bending  moment  expressions  obtained  using  the  bending 
moment  displacement  relations,  are  given  by  Equations  (44) 
and  (45).  Maximum  bending  moments  were  found  to  occur 
at  the  center  of  the  plate,  and  were  obtained  as  Equations 
(46)  and  (47). 

The  expression  for  maximum  deflection,  and  maximum 
bending  moments  in  x and  y directions  at  the  center  of  the 
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plate  were  found  to  vary  with  b/a,  the  plate  aspect  ratio. 
They  were  also  found  to  be  rapidly  convergent  series,  with 
convergence  to  the  exact  solution  obtained  for  the 
maximum  deflection  with  the  use  of  only  two  terms  of  the 
series  m = 1,  3,  and  convergence  to  exact  solution  obtained 
for  the  bending  moments  with  the  use  of  only  three  terms  of 
the  series  (m  = 1,  3,  5).  The  maximum  deflections,  and 
maximum  bending  moments  were  computed  for  varying 
aspect  ratios  (b/a)  ranging  from  1.0  to  4 using  three  terms  of 
the  series,  and  presented  in  Table  1.  The  table  shows  that 
the  coefficients  obtained  for  maximum  deflection  and 
maximum  bending  moments  were  exactly  identical  with  the 
solutions  obtained  using  the  Levy  single  Fourier  series 
method. 

VI.  CONCLUSIONS 

The  following  conclusions  can  be  made: 

i.  the  Kantorovich-Vlasov  method  can  be  successfully 
applied  to  the  solution  of  plate  flexure  problems  for 
simply  supported  ends  and  static  transverse  loads. 

ii.  the  solution  obtained  for  deflection  and  moments  is  a 
single  trigonometric  series  containing  hyperbolic 
functions,  which  may  not  be  so  easily  amenable  to 
mathematical  analysis. 

iii.  the  solutions  obtained  for  maximum  deflection  and 
maximum  bending  moment  are  single  trigonometrical 
series  that  have  rapidly  convergent  properties,  hence 
convergent  results  to  the  exact  solution  for  maximum 
deflection  and  maximum  bending  moment  are 
obtained  using  a few  terms  of  the  series. 
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